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T and J/^' spectroscopy using clover fermions in the presence of 
dynamical quarks 

S. Collins, R.G. Edwards, U.M. Heller, and J. Sloan^ 
^SCRI, Florida State University 

We calculate spin average splittings in the bb and cc systems using the clover action. We compare static and 
kinetic masses of heavy-heavy and heavy-light mesons and discuss the consistency of the results. 
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1. INTRODUCTION 

In this poster presentation, we discuss recent 
preliminary calculations of spectroscopy in the bb 
and cc systems, which were carried out using the 
tadpole improved Clover action and interpreted 
in the limit of heavy quarks using the Fermilab in- 
terpretation of this action [ |^ . In the first half, we 
discuss our results for the spin-averaged IP — IS 
and 25* — 15 splittings, which are commonly used 
to extract ajjg{Mz), and compare to those of 
other groups. In the second half, we investigate 
the relationship between the static and kinetic 
masses of mesons, Mi and M2 respectively, com- 
posed of various mixtures of 'heavy' and 'light' 
quarks (i.e. heavy-heavy, heavy-light, and light- 
light). The Fermilab interpretation of the Clover 
action with temporal and spatial k values equal 
(i.e. symmetric n) identifies the kinetic mass as 
the physical mass and differences between kinetic 
and static masses as a shift in the energy zero of 
the quarks. This interpretation requires 



dM{HH) + 6M{LL) = 26M{HL), 



(1) 



where SM = M2 - Mi, and HH, HL, and LL are 
different flavor combinations of heavy and light 
quarks. We find that this is obeyed, at this par- 
ticular value of n/ and /3, for M2(HL) < 1, but is 
strongly violated for M2{HL) > 1. The implica- 
tion of this is that, for example, there is no choice 
of Kb (for our ensemble) for which Af2(T) and 
M2{B) simultaneously agree with experiment^]. 



^This is due to the large dimensionless mass, not quench- 
ing effects on the running of the strong couphng 



2. THE SIMULATION AND FITTING 

The gauge configurations, smearing proce- 
dures, and fitting procedures are discussed in [ 
|j. For P states, an Z = 1 bound state spatial 
smearing function was used; one of the quarks was 
smeared by a gaussian multiplied hy x-\-iy, where 
X and y are the values of the x and y coordinates, 
and the other quark was a point source. For the 
"local" P-wave sink we used similar smearings 
which combined derivatives of a delta function. 
We used three different smearing radii each for 
bb and cc, and performed both matrix and vector 
fits with one, two, and three exponentials. All 
fits were done with tmax = 15, i.e. ommiting the 
central point of the propagator. 

Because onia (especially T) are so small, our 
physical box contained several relatively decor- 
related volumes. We took advantage of this by 
superposing the quark propagators from 8 spa- 
tial origins, using Z{Z) noise as described in [ ||. 
For P states, we further doubled our statistics by 
using t = 1& (in addition to i = 0) as an initial 
timeslice. This was vital in obtaining a reliable 
IP-IS splitting; the effective mass plots with the 
original statistics showed a "bump" near the cen- 
ter which led to ambiguous fits. 

For our determination of kinetic masses, we fit 
the ratio of a momentum one state and a zero 
momentum state to a single exponential, using a 
plateau of 8-12 in all cases. Looking at effective 
mass plots of the ratio, we saw that this corre- 
sponded to a very conservative value for imm, so 
our statistical errors should be much larger than 
the (systematic) fitting errors. 
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3. SPECTROSCOPY RESULTS 



4. KINETIC MASSES 



We present the results of a representative sam- 
ple of our fits in Table |^. The bb fits seemed better 
behaved than the cc fits; this can be seen by com- 
paring hb to he results. We see that corresponding 
T and r/f, fits are very similar; this is also true of 
the J ftp and rjc. We quote preliminary splittings 
m Table |; the error bars include our estimate 
of sytematic fitting errors. This leads to inverse 
lattice spacings of 2.27(13) and 2.l7{25)GeV from 
the bottomonium IP-lS* and 25-15' splittings, re- 
spectively, and 2.0(2) and 2.3{5)GeV from char- 
monium. For the respective dimensionless kinetic 
masses, aM2, of the lowest energy vector states 
we find 4.45(6) and 1.44(3) (corresponding to di- 
mensionful masses of 10.1(6)GeV^ for "T" and 
2.9(3)G'ey for "J/^f" when combined with the 
corresponding IP- 15 a~^.) Hyperfine splittings 
were estimated by examining effective mass plots 
of the ratio of vector to pseudo-scalar propaga- 
tors; rough values of Al{b)MeV for bottomonium 
and 80(20)Mey for charmonium were obtained. 

Table 2 



Preliminary 


sphttings 




hc{lP) - 


cc(15) 


0.23(2) 


J/tP{2S) - 


J/^(15) 


0.26(5) 


hb{lP) - 


56(15) 


0.199(11) 


T(25) - 


T(15) 


0.26(3) 



NRQCD bb and Wilson cc spectroscopy has 
been studied on this ensemble by other groups [ 
I, |._The NRQCD inverse lattice spacings from 
the bb IP-IS and 25-15 splittings are 2.44(7) 
and 2.37(10), respectively, roughly 10% higher 
than our values. This is consistent with experi- 
ence from quenched calculations, where NRQCD 
yields splittings about 10% smaller than corre- 
sponding clover calculations [ ||. Our cc IP- 
15 inverse lattice spacing is consistent with, but 
larger than that quoted for Wilson fermions in [ 
1^. Because the clover term mainly affects spin- 
dependent splittings, we expect agreement be- 
tween Wilson and clover calculations of the spin- 
averaged spectrum. 



QCD is a Lorentz invariant theory. Therefore 
the dispersion relation of physical states should 
obey the relation E{p^) = + A-P. In [ it 
was pointed out that, for non-relativistic parti- 
cles, a discretization which subtracts off a zero of 
energy for each (heavy) quark captures the essen- 
tial physics: 



(2) 



where the Ihs is typically called the static mass. 
Ml, and 6Mi is the shift in the zero of energy of 
the i'th quark in the state. The kinetic mass is 

defined as M2 — (2 '^^^^ ^ \p2=o)~^, and is equal 
to M for the dispersion relation (||). A useful 
definition of M2 given the splitting 5E = E{p'^) — 
E{0) is obtained by solving Eq. (||) for M; one 
obtains M2 = {p^ /25E) - 5E/2. 

The sum over energy shifts in Eq. (^) can 
be calculated in simulations; it is just the 
discrepancy between kinetic and static masses 
5M {state) = M2{state) — Mi{state). Because 
the energy shifts in Eq. (||) depend only on the 
quark content of the state, states with the same 
quark content should have the same 5M] the en- 
ergy shift associated with a B meson should be 
the average of the T and tt shifts, i.e. half the 
energy shift of the T. This leads us to the rela- 
tion 25M{HL) = 5M{HH) + 5M{LL) as a non- 
perturbative test of the consistency of the idea of 
attributing discrepancies between Mi and M2 to 
a shift in the zero of energy. 

The energy shifts 5Mi have been calculated for 
free quarks in [ pfl ; this leads to the relation 



2 ,pert (Ml 



e*-fi sinh(Mi) 



sinh(Mi) + l '-^^ 

In Fig. |l|, we have plotted our results for M2 (P5) 
vs. Mi{PS) for both mesons composed of degen- 
erate quarks (HH) and those of mixed mass(IIL). 
It is striking that both HH and HL mesons ap- 
pear to lie on a universal curve for this range 
of mass values. We also include the 'perturba- 
tive' curve for HH mesons M2 = 2M2.j,ert(Mi/2), 
which seems to describe the curve fairly well. In- 
dications of this behavior were first seen in [ |^ , 
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Table 1 

Selected bb and cc fit results. 



State 


type 


Smears 


Ncosh 




T 

max 


Q 


El 


E2 




VbinS) 


Vec 


a,b,c 


3 


5 


15 


.06 


2.4539(8) 


2.709(15) 


2.94(3) 




Vec 


a,b 


2 


10 


15 


.37 


2.4532(9) 


2.763(17) 






Vec 


b,c 


3 


2 


15 


.03 


2.4539(7) 


2.690( 7) 


3.19(6) 




Mat 


a,b 


3 


3 


15 


.11 


2.4543(7) 


2.70 (2 ) 


2.96(4) 




Mat 


a,b 


2 


8 


15 


.18 


2.4545(6) 


2.751( 7) 






Mat 


a,b 


2 


11 


15 


.72 


2.4533(6) 


2.701(15) 




T{nS) 


Vec 


a,b,c 


3 


5 


15 


.07 


2.4707( 9) 


2.714(15) 


2.95(3) 




Vec 


a,b 


3 


2 


15 


.29 


2.4691(13) 


2.76 (3 ) 






Mat 


a,b 


3 


3 


15 


.13 


2.4710( 8) 


2.71 (2 ) 


2.97(4) 


hb{nP) 


Vec 


a,b,c 


3 


1 


15 


.30 


2.666(7) 


2.87(2) 


3.083(19) 




Vec 


a,b,c 


2 


7 


15 


.41 


2.668(9) 


2.95(3) 






Mat 


a,b,c 


2 


9 


15 


.34 


2.669(8) 


2.94(4) 




r]c{nS) 


Vec 


a,b,c 


3 


6 


15 


.24 


1.2156(17) 


1.48(8) 


1.78(15) 




Vec 


a,b,c 


2 


9 


15 


.56 


1.2150(13) 


1.53(3) 






Vec 


a,b 


2 


8 


15 


.55 


1.2157(12) 


1.55(2) 






Mat 


a,b 


3 


3 


15 


.14 


1.2137(15) 


1.47(7) 


1.77(11) 




Mat 


a,b 


2 


9 


15 


.45 


1.2138(13) 


1.49(3) 




J/tp{nS) 


Vec 


a,b,c 


3 


4 


15 


.20 


1.2547(15) 


1.53(2) 


1.96(4) 




Vec 


a,b 


3 


4 


15 


.16 


1.2541(2 ) 


1.51(4) 


1.86(3) 


hc{nP) 


Vec 


a,b,c 


3 


2 


15 


.66 


1.481( 8) 


1.86(3) 


2.52(2) 




Vec 


a,b,c 


2 


7 


15 


.63 


1.459(18) 


1.79(7) 






Vec 


a,b 


2 


2 


15 


.52 


1.490( 5) 


1.863(10) 






Mat 


a,b 


2 


5 


15 


.27 


1.484( 9) 


1.79(4) 






Mat 


a,c 


2 


6 


15 


.14 


1.469( 8) 


1.84(2) 





slightly too heavy) . For the same choice of kj, we 
find M2{B) = ^A{'i)GeV - much hghter than the 
physical B meson! 

5. CONCLUSIONS 

Our results for spin-averaged splittings of the 
bb and cc systems showed no surprises. All gave 
inverse lattice spacings between 2.0 and 2.3GeV, 
with the smallest inverse lattice spacings com- 
ing from the softest splittings. This is consistent 
with a partial removal of quenching effects by the 
two flavors of dynamical fcrmions. The T inverse 
lattice spacings were roughly 10% smaller than 
in corresponding NRQCD calculations, consistent 
with experience in quenched calculations. 

A problem arises, however, when we exam- 
ine the relationship between static and kinetic 
masses. We find that, after tuning Kb so that 



albeit with large error bars. 

In Fig. 1^ we plot the relative inconsistency 

SAIjHL) - {SM{HH)+6M{LL))/2 ^^^ 
= AUHL) 

as a function of M2{HL) for many different 
heavy- light mesons. For M2{HL) < 1, we see 
that Eq. (|^) is consistent; for M2{HL) > 1, how- 
ever the relative discrepancy in energy shifts ap- 
pears to be growing linearly in M2{HL) for fixed 
light quark mass. Since the D meson should have 
a dimensionless mass of about .7 on this ensemble, 
the energy shift ansatz is consistent for charmed 
mesons. The B meson, however, should have a 
dimensionless mass of about 2.5; if we tune Kt, 
by matching M2(T) to experiment (as we did), 
then we expect M2{B) to be much too light. 
Recall that we found IvhiT) = 10.1(6)Gey for 
our choice of Kb (i.e. we used a bare quark mass 
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Figure 1. JVhiPS) vs Mi{PS) for HH and HL 
mesons. The perturbative curve for HH mesons 
is shown as a sohd hne. 
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Figure 2. I{H,L) vs M2{HL). Points with the 
same symbol have the same hghter quark mass. 

the kinetic mass of the T agrees with experiment, 
the kinetic mass of the B meson is roughly 20% 
too small. This will obviously cause difficulties 
with trying to calculate fs directly at the B me- 
son mass on a 2GeV lattice, especially since fs is 
strongly dependent upon Mb ■ This problem does 
not appear to be present in NRQCD; Eq. (||) has 
been checked for the Be system [ || . Another area 
where this inconsistency may cause problems is in 
low-/3 calculations; aMp w 1.5 when the inverse 
lattice spacing is 500MeV. With this in mind, 
the dependence of the inconsistency threshhold 
upon lattice spacing should be checked. 



Recently, another group has postulated a par- 
ticular form of the lattice dispersion relation for 
mesons [ ^, which leads them to similar conclu- 
sions to those presented here. We stress that we 
calculate both M2 and Mi directly from the sim- 
ulation, rather than relating them through a par- 
ticular ansatz for the dispersion relation. 

Finally, we would like to stress that our results 
do not prove that B physics cannot be done on a 
2GeV lattice. First, we know that the clover ac- 
tion does not include retardation effects correctly 
for very heavy quarks; it is possible that adding 
improvement terms to correct this would reduce 
the inconsistency. A more promising idea is to 
explicitly break the hypercubic symmetry of the 
clover action by using different k's in the temporal 
and spatial directions [ |^ ; the asymmetry could 
then be tuned to require that Mi agree with M2, 
obviating the need for energy shifts. 
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